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Abstract 

The continuous wavelet transform is adapted to account for sig- 
nal truncation through renormalization and by modifying the shape 
of the analyzing window. Comparison is made of the instant and in- 
tegrated wavelet power with previous algorithms. The edge adapted 
and renormalized admissible wavelet transforms are used to estimate 
the level of solar magnetic activity from the sunspot record. The solar 
activity is compared to Oerlemans' temperature reconstruction and 
to the Central England Temperature record. A correlation is seen for 
years between 1610 and 1990, followed by a strong deviation as the 
recently observed temperature increases. 

1 Introduction 

The continuous wavelet transform has become a popular tool in the analy- 
sis of data signals [U El El HI El El [7] , particularly among those investigat- 
ing variations in global climate [H El [TO] . Essentially just a special type of 
windowed Fourier transform, its primary drawback is the loss of transform 
response attributed to the finite length of a truncated signal, known as the 
cone-of-influence. Various algorithms have been proposed to correct for its 
effect [TT1 [T2l [T3| [H] ; however, none have proven entirely satisfactory. In this 
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article, we investigate the proper renormalization for a truncated wavelet as 
well as propose an algorithm which modifies the shape of the analyzing win- 
dow with respect to both the location of the wavelet and the length of the 
signal. These algorithms produce a reasonable estimate for the integrated 
wavelet power even as they approach the edges of the known data, verified 
by extensive evaluation of test signals. 

Armed with these transforms, we calculate the level of solar magnetic 
activity as given by the sunspot record using both the international [T5] and 
group [16] sunspot numbers, consistent with our previous estimate [H] . Our 
indicators of climate change are the global temperature reconstruction by [17] 
and the Central England Temperature record [18J. Plotting the temperature 
against the solar magnetic activity reveals a correlation which extends up 
until about 1990, after which a strong deviation from the linear trend is noted 
as the recently observed temperature increases despite a decreasing level of 
solar activity. The most recent year for which data is available reapproaches 
the trend. 

2 Normalization, Central Frequency, and Ad- 
missibility 

We begin by writing the complex Morlet wavelet ip® t (t') in its conventional 
normalization at scale s = l/f s = 2tt/uj s and offset t with unity sample rate 
At' = 1 as the product of a scale dependent constant C, a normalized window 
$, and a normalized wave using the parameter r\ = (t' — t)/s, 

< t (f) = CS$2 )t (Oe2 ft (f) = ir-^s-^e-^e^ , (1) 

where u\ ~ 2-k is the central frequency of the mother wavelet at unity scale 
and zero offset, V'ioCO = t 1 " ^ 4 e~' l 2 e iwxt . In practice, the window has a 
finite extent of — [sx <t'< [sx and the wavelet a length of N t > = 2 [sx + 
1, where the parameter x defines the resolution width beyond which the 
Gaussian is taken to be negligible; here we use x — 6- The mother wavelet is 
normalized to unit energy q {uj)\ 2 du) = 1 so that the Fourier transform 

of ip°(t'/s) is ip°(su) = \/27rs -01 o ( SUJ ) ■ We have found previously that the 
inclusion of additional factors in the scale dependent constant C s = a/2C!?/s 
produces a transform with some very desirable properties, 

^ S)t {t') = v / 2 7r- 1 /V 3/ V 2/ V^ . (2) 
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The factor a/2 comes from absorbing the transform response at negative 
scales, usually neglected, so that the mother wavelet now has a norm of 2, 
and the factor s~ 3 / 2 arises from symmetrization of the forward and inverse 
wavelet transforms of a signal y(t) = Re(A y e luJyt ) with duration N t by pulling 
over a factor 1/s usually left in the denominator of the integrating measure 
(cf. Equations (6) and (9) by [12]), 

CWT(m)= J2^At')y(t') , (3) 



ICWT(t) = Re 



(4) 



noting that = ip(—rf), so that the scaled wavelet now has a norm of 

2/s 2 . This normalization produces a power spectral density PSD(s, £) = 
|\/2CWT| 2 such that the integrated area of an isolated peak in the instant 
wavelet power IWP t (s) = PSD(s, t) returns the square of the amplitude A y 
of the signal element. Noting that the window has imparted the particle-like 
property of locality, we appeal to the photon whose energy is proportional 
to its frequency E u oc v = l/s u yet must extend for at least one cycle period 
to write its power as oc suggesting a physical basis for our choice of 

normalization. The factor \/2 applied to the amplitude of the CWT when 
forming the PSD we attribute to a root-mean-square effect, such that if 
we asked for the rms PSD, no factor would be necessary here but would 
then be introduced to compare the estimate for the rms amplitude to the 
peak amplitude A y . The margins of the PSD yield the mean wavelet power 
MWP(s) = Nf 1 ]T t PSDO,t) and the integrated instant power IIP(t) = 
J2 S PSD(s, t)As. The analysis of a test signal of duration 300 time units with 
signal elements of amplitudes 1 and v2 and periods 10 and 50 respectively is 
shown in Figure [H The cone-of-influence is defined as the e-folding time, and 
the more restrictive cone-of-admissibility indicates where wavelets are wholly 
contained by the extent of the signal. Apparent is the loss of transform 
response and reconstruction as the signal edge is approached, and the slight 
overestimation of the power is attributed to the response at large scales to 
the signal truncation. 

The central frequency given by |4] to unify the Fourier period and wavelet 
scale, Ai/si = 4tt/(u>i + a/2 + oof) = 1 yielding u>i = 2tt — 1/Att, is no 
longer appropriate for our normalization. Using a test signal of the same 
duration with elements of unit amplitude and periods 5, 15, and 50, we 
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consider transforms with central frequencies 2ti — l/4ir, 2tt, and 2ir + 1/47T 
and forward transform scalings of s -1 / 2 , s -1 , and s~ 3 / 2 appearing in the 
CWT. The left column in Figure [2] displays the instant wavelet power at 
the center of the transform N t /2 for a single central frequency (the others 
are shifted by a small amount), and the right column shows its gradient for 
all three central frequencies in the vicinity of the signal peak at period 15; 
similar graphs obtain for the other peaks, noting that the forward scaling of 
in Figures WLc) and (d) corresponds to that recently proposed by [19J. 
Only the transform with scaling s -3 / 2 produces peaks with an integrated 
area equal to the sum of the squares of the amplitudes, and we note that 
the locations of its peaks coincide with the signal periods for the central 
frequency of U\ = 2tt + 1/47T. 

In order to satisfy the admissibility condition ipi(0) = 0, one must use 
the zero mean formula so that the wavelet has no zero frequency response. 
Neglecting normalization and offset, the Fourier transform 

POO 

V>?M = / ^(t^e-^dt oc e ~^-" )2/2 (5) 



indicates the mother wavelet contains a DC component V'i(O) °c e^i/ 2 = 
d c ~ 10 -9 . From 0° = e luJlV one must subtract the continuous admissibility 
term, so that the zero mean wave becomes = O° — d c , defining the admis- 
sible wavelet transform AWT. With this subtraction, the Fourier transform 
of the admissible mother wavelet acquires a term which cancels the DC com- 
ponent 

/oo 
C^WOiWe-^'dt' , (6) 
-oo 

K _ e -KW)/2 ^ (7) 

and setting its derivative with respect to uj\ 

oc ( W - ^)e~(— ) 2 / 2 + u ie ~^ + ^ 2 (8) 

dui\ 

equal to to reveal the extrema yields the equation 

u 1 -w = cute'" 1 " , (9) 

where the effect of the zero mean formula d c is carried through on the right 
hand side. Its solution is given graphically in Figure [3] alongside u)\ = ui, and 
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its interpretation in the first and third quadrants as the wavelet response to 
positive and negative frequencies is straightforward, less so for the second 
and fourth where the second derivative indicates a minimal response. The 
solution crosses the u% axis at ±1, and the effect of the admissibility term is 
to induce a nonlinear response as the signal frequency approaches zero. 

Details of the edge adaption have been relegated to the Appendix. Our 
previous algorithm, identified as the corrected continuous wavelet transform 
CCWT, is found to overestimate the power carried by low frequency compo- 
nents. The newly proposed edge adapted EAWT and renormalized admissi- 
ble RAWT transforms perform much better at estimating the power carried 
by the signal. This performance comes at a cost to the reconstruction for 
the EAWT, as a trade-off between frequency response and power estimation 
has been made. 

3 Data Selection 

Turning now to the data chosen for this analysis, our indicator of solar mag- 
netic activity will be the integrated instant power of the yearly sunspot num- 
ber. As it is the underlying magnetic activity characterized by the Hale cycle 
of ~22 years in which we are interested, we must derectify the signal given by 
the sunspot number; doing so has the advantage of keeping the cycle shapes 
unaltered compared to simply analyzing the mean-subtracted rectified sig- 
nal. Previously, we followed [201 EI] in taking the square root of the sunspot 
number and alternating the sign as appropriate, but after a thorough eval- 
uation of the arguments presented by [221 [23], we now only alternate the 
sign to let any nonlinearity show up naturally as harmonics in the spectrum 
and subtract the mean before entering the wavelet analysis. We will con- 
sider both the sunspot group number R g by |16j . which covers the years 1610 
through 1995, as well as the international sunspot number Ri by the [15], 
which continues the Wolf index from 1700 up until the present day. While 
R g starts earlier and displays a more pronounced secular trend [21], Ri is 
known to display a tighter correlation with sunspot area and radio flux [25J ; 
their relative utility depends on whether one is more interested in the past 
or the future. The derectified signals are shown in Figure H^a). 

Our indicators of climate change are the glacial reconstruction of global 
temperature by [17] and the Central England Temperature record [18], "the 
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longest instrumental record of temperature in the world. "0 The glacial tem- 
perature reconstruction GTR spans the years 1600 through 1990 and com- 
bines data from around the globe, while the CET data begin in 1659 and 
continue to the present for a specific locality, both displayed in Figure IU(b) 
with their mean subtracted. Also shown is the 9-point running average of the 
CET, which generally oscillates about the GTR. The recent warming trend 
is apparent in both data sets, and the latest data indicate a downturn in 
observed temperature. Note that the GTR is defined as the "anomaly" and 
has had a constant subtracted by its originator. 

4 Analysis of the Data 

Beginning with the signal for Ri taken up until 2008, we display its EAWT, 
RAWT, and CCWT in Figure El The quality of reconstruction is not as good 
for the EAWT, as measured by A IWT the absolute value of the discrepancy 
with the signal, which is several factors larger than for the other transforms. 
The IIP for the CCWT generally agrees with that for the RAWT yet exceeds 
it for recent years, and the EAWT has a slightly different behavior for decades 
at either edge of the data. This total power estimate compares well with that 
found by [26] determined from just the Schwabe cycle. The spectra given by 
the MWP are in close agreement, except for a region around the 100 year 
scale where the response of the EAWT has dropped, with the third harmonic 
at ~7 years clearly present. The trough marking the extremely low frequency 
ELF region is not located precisely at the scale of the duration. Discarding 
the CCWT, the results for R g are shown in Figure [6j The IIPs of the RAWT 
and the EAWT again display some variation for recent decades. The MWP 
for R g does not display as prominent a third harmonic, and again the response 
of the EAWT is slightly diminished for scales between the maximum scales 
of admissibility and influence. No trough appears at the duration scale, 
indicating significant ELF contribution. We take the average of the IIPs for 
the EAWT and the RAWT integrated to the duration scale as our indicator 
of solar magnetic activity SMA for both Ri and R g ; the ELF contribution is 
generally less than 1 percent of the total. 

Finally we are in position to compare the level of solar magnetic activity 
with our indicators of climate change. To produce our comparison, we plot 
the time series of temperature in Figure b) against the time series of solar 

1 Met Office Hadley Centre for Climate Change 
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magnetic activity from Figures [5(e) and [6(c) and project along the time 
axis. These results are similar to our previous work [13], but we feel that 
here we have a better estimate of the SMA; the similarity indicates that 
the result is robust. In Figures [7(a) and (b) we plot the global temperature 
reconstruction GTR against the SMA for R g and Ri respectively. No relation 
is seen until the warming trend of the 20th century begins, then both the 
GTR and SMA increase to their highest levels. We learn from [TTJ that 
measurements for the GTR before 1900 becomes increasingly scarce, thus 
the lack of a visible relationship for the first three centuries may result from 
an inadequate supply of data. When the CET is plotted against the SMA 
in Figures [7J(c) and (d), we find a band emerges within which the points 
are scattered. A linear fit to the points is shown by the solid line, and the 
quality statistic Q = erfc(|rp| a/ N/2) determined by Pearson's r as well as 
the baseline temperature TJ, at SMA are displayed in Table [Hand compared 
with values obtained when applying the square root to the sunspot number 
as previously. Interestingly, the warming attributed to the increase in SMA 
by R g is .83 degrees, nearly the same as the total warming in the GTR 
of .78 degrees. The relation becomes more apparent when one uses the 9- 
point running average of the CET, Figures [7(e) and (f). Here we find the 
linear trend holds up until about 1990, after which the observed temperature 
increases without relation to the magnetic activity. We note that our results 
are in accord with other investigators [271 I2BT [TO] and contrast with those who 
have found negligible evidence of solar forcing on Earth's climate [29l [30l [3Tj . 



5 Discussion and Conclusions 

Many factors can influence global climate and local meteorology [32], [HI [331 El 
[331 [35], thus we should not expect a correlation with r P = ±1. The low value 
of Q indicates that the null hypothesis of no correlation should be rejected 
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in favor of the linear trend. Comparison to others' work is difficult as few 
authors make use of the derectified sunspot signal and most employ a cross 
coherence analysis of the spectra. The methodology here is quite different, 
eliciting a relationship between total solar magnetic activity and the observed 
temperature in central England. Naturally, investigation of the magnetic 
activity's correlation with other indicators of climate is forthcoming. The 
position of the most recent data point in Figure 0(f) indicates a possible end 
to the recent warming pattern displaying the most extreme excursion from 
the trend; only time can reveal whether our energy hungry civilization has 
averted a disaster beyond natural proportions. 

What influence has the derectification had on our results? From com- 
parison of the simple FFT spectrum of a sinusoidal signal with its rectified 
version in Figure [8], one finds that rectification turns a pure signal compo- 
nent into its even harmonic series and assigns most of the power to the DC 
component. Consequently, spectral analysis of the absolute sunspot number 
will contain overtones of the Hale cycle rather than its fundamental and at- 
tribute the majority of its influence to the constant component not addressed 
by the wavelet transform. The choice of where to apply the change in sign to 
the sunspot number generally spans only two or three years, thus introduces 
sporadic noise at the smallest scales which does not affect the overall result. 
Ideally, one wishes to work with a sunspot number which takes polarity into 
account, Rj = R + — where R± represent numbers of sunspots with op- 
posite magnetic orientation; however, such information is not available for 
the time span of the historical record. 

By adapting the continuous wavelet transform, one may estimate the in- 
stant power carried by a signal even as the edges are approached. Both renor- 
malization and window modification are found to be effective, where skewing 
the window improves the power estimation at the cost of frequency response 
and quality of reconstruction. Analyzing both the group and international 
sunspot numbers produces an estimate of the level of solar magnetic activity 
for each. When plotted against an ordinate of temperature, a relationship 
appears between the number of spots on the sun and the temperature ob- 
served in central England. This relation holds from the middle of the 17th 
century up until near the end of the 20th, after which the warming trend can 
no longer be explained by an increased level of solar magnetic activity. 
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6 Appendix 



6.1 Renormalization and Edge Adaptation 

Compensating for the effects of signal truncation is not straightforward. The 
CWT and AWT simply zero pad the signal beyond its known values, leading 
to the reduction in response at the edges. Previously [T3], we have sug- 
gested an algorithm to correct for the loss of response by renormalizing the 
CWT outside the cone-of-admissibility, and we describe it now in terms of 
a renormalization of the wavelet amplitude. For wavelets truncated by ei- 
ther edge of the signal, the window is shifted by an offset r relative to an 
unshifted window $o defining the time span t', and the corrected continu- 
ous wavelet transform CCWT is defined by C Sjt = C S ^JY S ^V, where 
Y $0 = s\/2tt for Y = Ye wh en written without subscript, and similarly for 
wavelets truncated on both sides. We have found lately that this algorithm 
may overestimate the power at low frequencies. 

As [12] point out, wavelet truncation affects the zero mean formula, such 
that the continuous admissibility term is no longer adequate to remove the 
DC component. The adaptive admissible wavelet transform AAWT is defined 
by 6 = 6° - d s>tl where d s ,t = (0°) = Y © 0( &r/ Y ®t is the weighted mean 
of the remaining wave 0°. Addressing now the amplitude of the truncated 
wavelet ip T = C S $ T 0, with C S)t = C s ^j2/ s 2 \ip T \ 2 we define the renormalized 
adaptive wavelet transform RAWT. These transformations are very similar 
to those encountered when building operators in lattice gauge theory [36] 
and in other applications. 

Finally, we consider an algorithm which modifies the shape of the an- 
alyzing window $ with respect to the location and length of the signal. 
The length of a truncated wavelet ip T is defined to be the lesser of the raw 
wavelet length or the signal length, N T = min(A r i /, N t ). The offset r(t) 
is determined from either the center of the signal or the location of the 
cone-of-admissibility, and the algorithm to keep everything aligned gets a 
bit complicated: for t' = max(0, [sx — \_N t /2) and t'(k) = [— [sx, [sx] in- 
dexed by k from 1 to N?, if r < then t' -> t' [1, min(A/" t /, iV t )] + r', else 
t' — > £'[max(l, Nf — N t + 1), Nf] — r 1 . The end result is simply to truncate 
either edge of the wavelet as necessary. Then, for 7 = min(|_sx, \_Nt/2), the 
edge adapted window is ln$f t (t') = —[777/(7 ± r)] 2 /2 as t' < or > r, respec- 
tively, where $f t (r) = 1, which scales the t' axis of the window to either side 
of r. These wavelets, renormalized as above, define the edge adapted wavelet 
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Table 2: Summary of the wavelet transforms 
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transform EAWT. In Figure 0(a) we show the adapted windows for a scale 10 
wavelet centered at the labelled ticks, and in Figures [3(b)- (d) we display the 
real part of that wavelet for the EAWT and the RAWT for several offsets. 

6.2 Comparison of the Adapted Transforms 

We next embark upon a lengthy exploration of the relative responses of the 
various transforms. The properties of our transform menagerie are summa- 
rized in TableEl Let us first consider the responses of the EAWT, the AAWT, 
and the CCWT to a succession of test signals yu of duration N t = 300 with 
squared amplitudes A\ = [1, 2, 3, 4] and periods Pk = [5, 30, 150, 300]. In 
Figure [10] we display the IIP for each test signal. Generally, the CCWT is 
found to overcompensate for the wavelet truncation while the EAWT dis- 
plays some oscillation in the IIP near the edge. A slight overestimation of 
the power is attributed to the signal truncation. The IIP for the AAWT 
drops significantly in response to the wavelet truncation, especially for signal 
periods wholly beyond the cone-of-inffuence, while the IIP for the EAWT re- 
mains fairly close to the signal power for the duration. In Figure [Til we show 
the IWP at the midpoint of the transform for each of the test signals. Here, 
the peaks differentiate as the signal period approaches the signal length. We 
notice in Figure flTT d) that the peaks are not located precisely at the signal 
period. Intrigued, we next consider a series of test signals of unit amplitude 
with periods Pk = [350, 400, 500, 600] exceeding the duration, Figure [T2l As 
the signal period increases, the magnitude of the response diminishes, and 
the location of the peaks seems bound from above by the scale of the du- 
ration Nt- Additional peaks appear at harmonics of the response for these 
extremely low frequency ELF signals. Such behavior is not unexpected in 



10 



light of Figure [3J where one sees that the wavelet response to ELF signals 
has a finite bound in scale. 

Finally, we consider the responses of the EAWT, the RAWT, and the 
CCWT to a test signal of the same duration comprised of the sum y = 
Ylt=i V k OI " the first three signals above, with their PSD displayed in Fig- 
ures 1131(a)- (c). The CCWT produces a deeper trough between the peaks 
in the MWP of Figure 1131(d) and a generally stronger response, while the 
EAWT returns a deeper trough at the scale of the duration. In Figures [TW a) 
and (b) we display the reconstructions for the newly proposed transforms. 
The loss of reconstruction for the IEAWT is of a different character than for 
the others, and we wonder if an alternate reconstruction basis would be more 
effective. The best reconstruction is provided by the IRAWT, which still dis- 
plays some edge alterations. The IIPs of Figure [IHfc) have been evaluated 
with a cutoff at N t , while those of (d) include the contribution from the ELF 
region. The effect of the renormalization on the RAWT is clearly seen, as 
now its IIP agrees with that of the EAWT at the midpoint N t /2, while the 
CCWT overestimates the signal power. The ELF contribution is apparent as 
an increase in the IIP above the instant signal power and is attributable to 
the effect of signal truncation. We conclude that the EAWT and the RAWT 
each provide a reasonable estimate of the instant signal power. 
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Figure 1 : CWT analysis of a signal with elements of amplitude 1 and y/2 and 
period 10 and 50. (a) The power spectral density with the cone-of-influence 
(dashed) and cone-of-admissibility (dotted) overlaid. An unlabelled tick on 
the scale axis appears at the signal duration N t = 300. (b) Mean wavelet 
power displaying the signal peaks as well as the large scale response to the 
signal truncation. A trough appears at the scale of the signal length, (c) 
The integrated instant power decreases to half its expected value towards 
the signal edge, (d) The reconstruction is diminished compared to the signal 
(dots) near the edge 
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Figure 2: CWT analysis of a signal with elements of unit amplitude and pe- 
riods 5, 15, and 50 for central frequencies 2tc — 1/Att (dash-dot), 2n (dashed), 
and 27T + 1/Att (solid) and forward transform scalings of s~ 1//2 , s" 1 , and s~ 3 ^ 2 
labelled by the vertical axis. The left column is the instant wavelet power 
at N t /2 and the right is its gradient Ifor all three central frequencies in the 
vicinity of the signal peak at period 15. The area under each peak in (e) is 
unity 
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Figure 3: The solution (solid) of the equation uie^ 1 ^ = uj\ — u derived by 
setting dipiiu) / duJi = crosses the u)\ axis at ±1. The inadmissible solution 
(dashed) is linear 
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Figure 4: (a) Derectified signal for the group sunspot number R g (dashed) 
and the international sunspot number Ri (solid), (b) Mean subtracted glacial 
temperature reconstruction GTR (dashed) along with the Central England 
Temperature CET (dots) and 9-year smoothed CET (solid) 
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Figure 5: International sunspot number Ri power spectral density (a)-(c), 
quality of reconstruction (d), integrated instant power (e), and mean wavelet 
power (f) for the EAWT (dashed), the RAWT (solid), and the CCWT (dash- 
dot). A IWT is the absolute value of l2ie discrepancy between the reconstruc- 
tion and the signal 
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Figure 6: Group sunspot number R g power spectral density (a)-(b), in- 
tegrated instant power (c), and mean wavelet power (d) for the EAWT 
(dashed), the RAWT (solid) 
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Figure 7: Scatter-plot of the glacial temperature reconstruction (a)-(b) and 
the Central England Temperature retgg'd (c)-(f) versus solar magnetic activ- 
ity for the group number (left column) and the international number (right 
column). Shown is their projection along the time axis. The dotted line 
on the CET plots displays a linear fit, and CET 9 uses the 9 point running 
average to reduce the scatter. The data for R g end in 1995, and the point 
for year 2008 is marked with an x in (f) 
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Figure 8: Comparison of power spectrum of a signal (x) with its rectified 
version (+) both with and without subtraction of the mean. The power at 
frequency equals 0.405 
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Figure 9: (a) Edge adapted windows <& E for a scale 10 wavelet centered at 
the labelled ticks. The raw wavelet length is 121 units, (b)-(d) Real part 
of the scale 10 wavelets for the EAWT (dashed) and the RAWT (solid) for 
several offsets 
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Figure 10: Integrated instant power for the EAWT (dashed), the AAWT 
(solid), and the CCWT (dash-dot) of a test signal y P = Re(e it2w/P ) with 
period P 
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Figure 11: Instant wavelet power at N t /2 for test signals with period P 
indicated by the vertical dotted line for the EAWT (dashed), the AAWT 
(solid), and the CCWT (dash-dot). Note that the y-axes are in milliunits. 
In succession, the integrated power for the EAWT is 1.01, 2.04, 3.08, and 
3.84 
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Figure 12: Instant wavelet power at the midpoint of the duration for test 
signals with period P for the EAWT (dashed), the AAWT (solid), and the 
CCWT (dash-dot). Note that the y-axes are in milliunits. In succession, the 
integrated power for the EAWT is 0.97, 0.78, 0.34, and 0.15 
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Figure 13: Power spectral density (a)-(c) and mean wavelet power (d) of the 
test signal y = y 5 + y 30 + yi 50 for the EAWT (dashed), the RAWT (solid), 
and the CCWT (dash-dot) 
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Figure 14: Reconstruction from the IEAWT (a) and the IRAWT (b) and 
integrated instant power (c)-(d) for the EAWT (dashed), the RAWT (solid), 
and the CCWT (dash-dot). The IIP of (c) has been cutoff at the scale 
N t = 300, while that of (d) includes the ELF contribution 
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